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Abstract 



We present a calculation scheme for the two-loop vacuum polarization cor- 
rection of order a 2 to the Lamb shift of hydrogen-like high-Z atoms. The 
interaction with the external Coulomb field is taken into account to all or- 
ders in (Zct). By means of a modified potential approach the problem is 
reduced to the evaluation of effective one-loop vacuum polarization poten- 
tials. An expression for the energy shift is deduced within the framework 
of partial wave decomposition performing appropriate subtractions. Exact 
results for the two-loop vacuum polarization contribution to the Lamb shift 
of K- and L-shell electron states in hydrogen-like Lead and Uranium are 
presented. 



PACS-numbers: 31.10.+Z, 31.30.-i, 31.30.Jv 



I. INTRODUCTION 



Recent experimental progress in the spectroscopy of highly charged heavy ions JTJ-Q de- 
mands theoretical predictions for the Lamb shift, which should include the complete set of 
QED radiative corrections of order a 2 but accounting for all orders (Za) in the interaction 
with the strong external Coulomb field. For \ow-Z elements, a potential expansion with 
respect to powers in (Za) is legitimate and all a 2 -corrections have recently been calcu- 
lated up to the order of a 2 (Za) 5 |j]J|. However, for systems under consideration in recent 
Lamb-shift measurements such as Gold or Uranium, a value for the effective coupling 
Za > 0.5 already indicates that Za-expansion in the regime of large Z becomes inade- 
quate. Instead, exact electron propagators and wave functions in the external Coulomb 
field of extended nuclei have to be used in calculations of all second-order diagrams (Fig. 
|]). Meanwhile, calculational approaches are available for most of these QED effects (see 
the work of Persson et al. |7j and cited references). However, the calculations are not yet 
complete: 1. Exact numerical evaluations of the complete set of two-photon self-energy 
contributions remains as a major challenge, although recent progress has been made in de- 
riving renormalized expressions for the resulting energy shifts [12| . The lack of numerical 
results for these contributions represents a major uncertainy in theoretical predictions for 
the Lamb shift aiming for a relative precision of 10~ 6 for the total electron binding energy. 
2. Exact evaluation schemes, which also treat the loops involved to all orders in (Za), 
have been developed only for the combined self energy - vacuum polarization correc- 
tions SEVPabc and the two-loop ladder vacuum polarization diagram VPVPa (see Fig. 
[I]). Values for the VPVPa correction for large Z-numbers have been tabulated recently 
||. The self energy-vacuum polarization S(VP)E (Fig. |l|) is calculated by employing the 
Uehling-approximation for the loop J7|. Calculation of the higher-order contribution to 
this diagram are presently in progress. 3. Until now, the two-loop diagram VPVPb as 
well as the self energy corrected one-loop vacuum polarization contribution VPVPc (Fig. 
[1]) where calculated only to lowest order in (Za) utilizing the Kallen-Sabry polarization 
function 0,|,[lO|]. 

Aiming towards the completion of the exact evaluation of all QED-radiative corrections 
of order a 2 , we shall present a calculation scheme which allows the higher-order (Za)- 
contribution of the two-loop vacuum polarization correction VPVPb to be determined. 
Recognizing that the diagram we wish to calculate is part of the complete " dressed" one- 
loop vacuum polarization allows us to reduce the problem to the evaluation of effective 



one-loop corrections where the renormalization procedure is known [|13| .|T5 



Section II contains a general discussion of "dressed" electron propagators in arbitrary, 
classical external fields and of the induced vacuum polarization. In section III we specify 
the formulae to the situation in hydrogen-like ions. The subtraction scheme for deducing 
the two-loop correction from various effective one-loop contributions will be introduced. 
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Section IV briefly reviews the renormalization procedure. An expression for the renor- 
malized energy shift due to the higher-order (in (Za)) part is derived. This consists of 
two terms which require a different numerical treatment. In order to compare the effect 
due to higher orders in (Za) we also present results for the two-loop correction in Uehling 
approximation, which will be derived in section V. In section VI we will calculate the 
effect of the higher-order contribution to the lS-Lamb shift in hydrogen-like Lead and 
Uranium. 

Throughout this paper, units will be used where H = m = 1 and e 2 = a. 

II. DRESSED ELECTRON PROPAGATORS AND ONE-LOOP VACUUM 

POLARIZATION 

We start with a brief discussion of the concept of " dressed" electron propagators and of 
the corresponding "dressed" one- loop vacuum polarization (VP). Some general formulae 
will be derived, which will be employed in the next section. 

We shall adopt the term "dressed" electron (positron) line ip for an electron (positron) 
moving in an arbitrary, external electromagnetic field A^. The wave function that ac- 
counts for the interaction with this external field is a solution of the Dirac equation: 

[i? x - e4(x) - m] *P(x) = . (1) 

Choosing free electron lines as a reference, we shall refer to the state ip as "A-dressed" 
electron line. Similarly, an electron interacting with the external Coulomb-potential 
V c = cAq generated by the (bare) nuclear charge density distribution may be called 
a "Coulomb-dressed" electron. In general situations it is appropriate to divide the total 
external field into two parts A c ^ and A^, i.e.: 

A,(x) = A;(x) + A,(x) , (2) 

where e.g. the second term A^ may be considered a perturbation. Accordingly, considering 
electron states in the external field A c ^ as unperturbed states, we may then call these 
states ip "A-dressed". The propagators Sp(x,x') and Sp(x, x') describing electrons in the 
external field A^ respectively A^ are defined by 

[\$ x - e4(x) - m] Sp(x, x') 
[i@ x -e4 c (x)-m]S F (x,x') 

Note that they also satisfy the equations 

S£(x,x') [i^, + e4(x') + m] 
Sp(x, x) [\$ x , + e4 c (x') + m] 



= 5(x — x') , 

= 5(x - x') . (3) 



= — 5(x — x') , 

= -5(x - x') , (4) 
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where the (adjoint) Dirac-operators are acting to the left. In analogy to the Dyson 
equation defining dressed propagators in terms of improper self-energy insertions etc. we 
postulate an equation of the form 

Sp(x,x') = Sp(x,x') + J d 4 xi d 4 x 2 Sp(x, xi)JC(xx, x 2 )Sp(x 2 , x') , (5) 

with a kernel /C. With the aid of equations (§) and @) we can solve for the kernel: 

K(x, x') = e4(x) 5{x - x) + e4(x) S$(x, x')e4(x') . (6) 

Insertion of Eq. (|6|) into (^) leads to an equation for the "A-dressed" propagator Sp 
(taking Sp as the reference) which can be solved iteratively. The exact propagator Sp 
appears as the sum of the unperturbed external field propagator Sp, a part describing 
a single-scattering with the additional external field and a higher-order part which 
accounts for multiple-scattering contributions. Given a representation for Sp a formal 
expression for the corresponding one-loop vacuum-polarization current induced by the 
total external field can be derived. For later purposes we will already specialize to the 
case of static external fields. Since the propagators are homogeneous in time, one obtains: 

= ie/^ l [ Tr[ 1 »S F (r,r,E)} + J d 3 n Tr [f£%(r, n, E) e4{r x ) S F {r u r, E) 
+ J d 3 n d 3 r 2 Tr [^S F (r, n, E) e^(n) «S#(n, r 2 , E) e4(r 2 ) S c F (r 2 , r, E)] } . (7) 

Putting aside questions about renormalization for a moment, this formally exact equation 
implies that the one-loop vacuum polarization induced by the total field A^ is given 
as a sum of three terms: a part induced by the external field A e , a single-interaction 
contribution and a third part taking into account multiple interactions with the additional 
external field A^. 

The representation of the propagator Sp and of the induced vacuum polarization J Ap " 
derived above are not unique. The reason for this is provided by the fact that the decom- 
position of the total external field (0) is completely arbitrary. In particular, we could have 
chosen the free-field configuration as unperturbed reference. Consequently, the external 
field propagator Sp in the defining equation @ has to be replaced by the free Feynman 
propagator Sp leading to a similar kernel (|B|) which will contain the total external field 
Afj,. In this case Eq. (|7]) takes the form 



J 



Am/ 



le 



f l/d 3 nTr 
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"S£(r-n,£) e4( ri )Sp( ri -r,E) 



+ / d 3 n d 3 r 2 Tr \YS F {r - n, E) e4{r x ) S A {r u r 2 , E) e4{r 2 ) S F {r 2 - r, E) 
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The free closed-loop contribution vanishes in accordance with the Furry-theorem. The 
induced vacuum polarization itself gives rise to a modification A^ of the total external 
field: 

A,(r) = J dV D, v (r - r', 0) J K \r>) . (9) 
The free photon propagator is given by (in Feynman gauge) 

D^r - r', 0) = 9lv D(r - r' , 0) = - 9iat [ e -Mr-r>) 1 

1 r°° 

= g, u - AkY,^Y lm {r)Yf m {r')u{kr)u{kr') . (10) 
vr Jo fz 



e,m 



The last line of the expression above specifies the partial wave decomposition of the photon 
propagator. A spherically symmetric, static external potential V = eA induces only a 
static vacuum polarization charge density J A0 which will also be spherically symmetric. 
It gives rise to an effective one-loop VP-potential: 

V A (r) = ^-Jdkj (kr)Jdr'r' 2 jo(kr')J AO (r') . (11) 



III. SUBTRACTION SCHEME 

The general considerations of the previous section may have already anticipated how 
we are going to deduce the two-loop vacuum polarization correction. The basic idea is 
to derive this contribution from an effective one-loop vacuum polarization V A , which is 
dressed with the renormalized first-order vacuum polarization potential induced by 
the external Coulomb field V c of the nucleus. At first we need to specialize Eq. ([?[) to (|TTD 
to the situation of bound-state QED. In the presence of a static, spherically symmetric 
nuclear charge density, we specify the total external field @ as the sum of the (bare) 
external Coulomb potential V c and the renormalized, first-order vacuum polarization 
potential : 

V(r) = eA (r) = V c (r) + v£(r) . (12) 
The one-loop potential is obtained from 

V r y n p (r) = elo(r) = ia J dV D(r - r>, 0) U ^ Tr [ 7 °^(r\ r' , E)] \ (13) 

^ / ren 

after renormalization. Sp denotes the electron propagator in the external Coulomb field. 



In view of Eqs. (0) and (JlTj) the VP-dressed one-loop vacuum polarization potential 
formally reads 
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V(r) 



dE 



dVD(r-r',fl) / — Tr j°Sp(r', r', E) 



ia / dVD(r-r',fl) 



2tt 
27 



Tr 



+ /dSTr 7 °^(r',r 1 , J E;) 7 Vr(r 1 )^(r 1 ,r',E) 



Ot/VP/ 



+ / d 3 nd 3 r 2 Tr 7 °S F V, r lf £) 7 lO 1 ) S£(n, r 2 , i^Cf (r 2 ) S£(r a , r', £) 



,Ot/VP 



Ot/VP/ 



(14) 



A graphical representation of the VP-dressed one-loop potential is depicted in Fig. |^. 
The two-loop vacuum polarization correction we are interested in appears as the second 
term of the decomposition (14) together with the ordinary one- loop vacuum polarization 
(first term). The third term accounts for all multiple interaction terms and is at least of 
order a 4 (Za) 3 , since the additional potential contributes with a leading order aiZa) 
(Uehling-term). Thus, it is legitimate to neglect this higher-order part. Consequently, we 
identify as the renormalized two-loop potential 



VPVPb/ 
rcn \ 



ia 



d 3 r'D{r - r',0) 



dE 
27 



Tr 



dE m 
Tr 

2tt 



j°S^(r',r',E) 
j S^r',r',E) 



(15) 



after appropriate renormalizations have been applied to the right-hand side of the above 
equation. Note, however, that the formal expression for the effective one-loop vacuum 
polarization density (first term in the curly brackets) already contains the renormalized 
one- loop potential V^f , assuming that the renormalization of the exterior VP-loops can 
be performed separately. This issue will be adressed in the next section. We are now in 
the position to write down the energy shift of a bound electron state 4>a due to the exact 
two-loop correction: 



AjE VPVPb 



(16) 



In order to deduce exclusively the contribution Ai?J PVPb (h.o.) arising from higher orders 
in (Za), one has to subtract the corresponding two-loop contribution evaluated in the 
Uehling-approximation. 



IV. RENORMALIZATION 



Taking the above into consideration, the problem of renormalization of the exact two-loop 
potential 



U 



VPVPb/ 



dE 



2tt 



r) = ia I dV%-r',0) / d 3 r x / — Tr 7 °S'p (r', r 1 , E) 7° Sp(r 1 , r', E) 



,0 cC, 



x / d 3 r 2 D( ri -r 2 ,0) 



dE' 



Tr 



7 °^(r 2 ,r 2 ,£') 



(17) 
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reduces to the problem of renormalization of the one-loop vacuum polarization. This 
is sugguested by Eqs. flU] ) and flT5|) . It relies on the fact, that the external VP-loop is 
properly taken into account in terms of a VP-dressed effective one-loop, where the external 
VP-loop may be replaced by the renormalized first-order potential V^f . This is supported 
by the notion that the electron (positron) experiences the effective nuclear charge modified 
by the induced vacuum polarization cloud. Procedures for renormalization of one-loop 



potentials involved in Eq. (15) are well known from the evaluation of the energy shift due 
to the first-order vacuum polarization in external Coulomb fields. We adopt the partial 
wave decomposition approach together with the subtraction scheme developed in [ p~3| , |T5]j 
and apply it to Eq. ([15]). In order to perform similar steps as in the case of the ordinary 
Coulomb-dressed VP-loop, we employ the equivalent representation according to Eq. (||) 

V{r) = ia J d 3 r'D(r-r',0) J ^ [J d 3 n Tr [-f°S®(r' - r u E) 7 ° V{r x ) S£(n - r', E) 
+ / d 3 n d 3 r 2 Tr \j°S%(r' - r u E) 7 ° V{r x ) $*(r u r a , E)f V(r 2 ) 5°(r 2 - r', E) 



where V denotes the total external potential Eq. ([121) . Fig. ||] shows the diagrammatic 
representation of Eq. (|T^). Further evaluation requires a partial wave decomposion of all 
the propagators 





r',E) 


= E 






n/!/i 




r',E) 


= E 








S° F (r- 


r',E) 


= E 









Vw( r )Vw 


(r') 


E £ nK ^{\ 


-ir)) 




r') 


E — E UjK (l — 


\rj) 




>') 


E — s PiK (l — 


irj) 



(19) 



together with Eq. (|10"D for the photon propagator. According to the renormalization 
prescription developed in []15] the renormalized one-loop potential V ren is obtained as 
the sum of the finite "Wichmann-Kroll"-type contribution V^2 and of the renormalized 
"Uehling"-type potential Vp\. We are lead to the following expressions: 



a f ' Kmax f 

VM r ) = / dk] {kr) \ E sig n (^n, K ) (Vwl jo(fcr') \^ nK ) 

!i |«|=1 1 " 



_4^^ (<P P k\ 3 o(kr') Wp'k){v p 'k\ V \<P P k) 



£p,K £p',K 



Vn{r) 



dr' Amr' ep (r') f(r,r') 



(20) 



together with the radial kernel 

2a r°° 



e- 2r Z sinh(2r'fl 

Br — r ) — 

r£ 2r'£ 
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„, , . e~ 2r 't sinh(2rO 
+6(r' - r) — - 



a 



3nr 



r'£ 2r£ 
' (jt) [X2(2|r-r'|)- X 2(2(r + r'))] , 

(21) 




p v denotes the sum of the (bare) nuclear charge density p nuc and of the renormalized first- 
order vacuum polarization charge density induced by the Coulomb field of the nucleus. 
Eqs. (E^) and (^TJ) are analogous to expressions one has to deal with when deriving the 
renormalized, exact one-loop potential (0). Performing now the subtraction as implied 
by Eq. (|T5|), we identify the renormalized two-loop vacuum polarization potential: 

Wl PVPb (r) = Wj PVPb (r) + Wj 2 PVPb (r) , (22) 

U£ YP \r)=f o dr'4vrr' 2 [ep^(r')] f(r,r>) , (23) 

U^ Ph (r)~--j dkj (kr) £ £ sign(5„, K ) (^ nK \ jo(Ar') |Vw> 



k|=1 

(wl jo(fcf') |yp'/c)(y P '/c| ^ lw) 



-4EE 

E sign(£„ >K ) (0„ K | jo(fcr') |0 nK ) - 4^Z! 



p p' £p,K £ P',K 

(wl jo(kr') \<Pp> K ){<p p > K \ V c \<p pK ) 
p p' ^p> k %>') K 



(24) 



The term W^f VPb denotes the finite (renormalized) Uehling-type potential generated by 
the Uehling part of the renormalized first-order vacuum polarization charge density. The 
term Z^PVPb 

summarizes the finite part of the Wichmann-Kroll-type potential. Since 
the sum over partial waves terminates at some finite |ft m ax| each contribution in square 
brackets becomes well defined. The procedure described in Eqs. ( f22|) to ( j2"4"|) for deducing 
the renormalized exact two-loop vacuum polarization potential is depicted in Fig. |j. 
For numerical evaluations, however, the representation fl2~4] ) of the potential W^ PVPb is 
ruther cumbersome since it involves terms given as multiple summations over the free 
Dirac spectrum. A representation which is more convenient for numerical calculations 
can be introduced as follows: At first, we combine the one-potential terms in Eq. fl2~4] ) 
to a one-potential term involving the renormalized one-loop vacuum potential V^ p only. 
Secondly, the resulting one-potential term is replaced by the completely VP-dressed (free) 
one-loop potential assuming that the effects due to multiple-interaction contributions, 
which are at least of order a(a(Za)) 3 , are neglegible. Thus we obtain the alternative 
representation: 
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u 



VPVPtw 
T2 \ 



/ dkj (kr) \Y1 sign(£„ A ) (i/; nK \ i (kr') \ip nK ) 

71 J \k\=1 I n 

sign(£„ ift ) ((p nK \ jo(kr') \<f> nit ) - ^ sign(£ n Ho{kr') \(p nK ) 



(25) 



where the states (p nK are solutions of the Dirac equation with the external potential V^f 
and corresponding energy eigenvalues e n>K . 

Having derived the renormalized potential W r ^ VPb , the corresponding energy shift of a 
bound electron state 4>a in the external Coulomb field can be evaluated according to (|16D . 

V. TWO-LOOP CONTRIBUTION IN UEHLING APPROXIMATION 



Being interested primarily in the contribution to the energy shift due to higher orders in 
(Za), we should subtract the Uehling-in-Uehling contribution, where the exact vacuum 
polarization loops are replaced by free fermion loops. Taking a uniform sphere model for 
the nuclear charge distribution, i.e. p nuc = (3Ze/4iiRl) 6(i?o ~ r ), the latter reads: 



7/VPVPb/ \ 



dr' 4nr 



ep r U e e n h (r') 



f(r, r') 



ep 



Ueh/ 



1 \ 1 



x le(Ro-r) 1 + 



sinh(2r£) e~ m ^ 



-e(r-Rn 



cosh(2i? 



sinh(2i? 



-2r£ 



The expression for p^ h may be cast into the more familiar form: 



ep 



Ueh/ 



— % ) 2- {sign^o - r) Xl (2\Ro - r\) - Xl (2(i? + r)) 



47ri?o / 



2Ro 



lX2(2\R -r\)- X 2(2(R + r))} 



(26) 



(27) 



This vacuum polarization charge density is plotted in Fig. |6]. It is easily verified that the 
induced Uehling density (p7|) possesses a logarithmic singularity at the nuclear radius, 
which originates from the first term in the curly brackets. Integrating the Uehling charge 
density over the interior of the nucleus we obtain the finite induced vacuum charge 

i-Ro 



eQ [nt (R ) = / dr'Airr' 2 
Jo 

Za a 



2Rq 71 



(Mm (at?\j- *3(4^o) , X 4 (4i? ) - X4(0) 
X2(0) + X2(4it ) H ^ h 



dr' 47rr 
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Ro 



R 

-eQ ext (Ro) 



(2R ) 



(2f 
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which is exactly cancelled by the total induced charge in the exterior region Q ext . 

The corresponding energy correction to the binding energy - we refer to it as Uehling-in- 

Uehling correction - which is part of the Kallen-Sabry correction |TT1 reads: 



AEY^l, Ueh) = (4> A \ ZYjPVPb ^ . (29) 
This will be calculated separately. 



VI. EVALUATION 

We now turn to the evaluation of the higher-oder (Za)-contribution of the two-loop 
vacuum polarization to the Lamb shift of strongly bound electrons 

A£7 VPb (h.o.) = (<f> A \ (^PVPb + w vPVPb) ^ m 

where the Uehling-in-Uehling part is subtracted. Thus the effect of all higher orders 
in (Za) in the interaction with the external Coulomb potential contributing to the exact 
two-loop correction proceeds in two separate steps. Accordingly, we only need to calculate 
the renormalized Uehling-potential generated exclusively by the induced Wichmann-Kroll 
charge density 

/■OO r -, 

Ujl^{r)= dr'^[epZ(r')]f(ry) , (31) 
Jo 1 J 

and the renormalized Wichmann-Kroll-type Potential W^ PVPb according to the subtraction 
scheme (E4|). 



As a first step towards the evaluation of the energy correction fl30|) we consider the con- 
tribution 

AEj PVPb (^l,WK) = (0 A |^Pb| 0j4) (32) 

This correction to the Lamb shift of the bound state (f>A is related to the change of the 
Uehling potential Eq. ( |5T| ) arising from the Wichmann-Kroll part of the induced vacuum 
polarization. At first we wish to derive an estimate for this correction to the Lamb shift 
of the ground state in hydrogen-like Lead and Uranium. 

The Wichmann-Kroll charge density is calculated based on the partial wave decom- 
position of the Coulomb propagator and of the free propagator as developed in Ref. |H 
It is obtained from 



n Jo 27r m=i Ui 

+ / dr'r' 2 V c {r') £ F^'l 
»,i=i 



r, r', iu) 



(33) 
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where the summation over k is terminated at some maximum value |ft m ax|- The quantities 
G™ and F l K J denote the partial wave decompositions of the free and the bound propagators. 
Apart from the long-range tail where the Wichmann-Kroll vacuum polarization charge 
density is positive, a strongly pronounced maximum of negative charge density occurs 
in the vincinity of the nuclear surface (see Ref. [13|]). Integrating over this r-range one 
obtains a total induced negative charge 



Q WK = eZ WK = ^ dr4ff/ ,WK (r) (34) 

Jo 

Since the Wichmann-kroll density is strongly localized near the nuclear surface it acts 
almost like an additional negatively charged spherical shell surrounding the nucleus. A 
strongly bound electron experiences the reduced nuclear charge. This suggests replacing 
the corresponding potential U^Ywk by the Uehling potential generated by a spherical 
shell carrying the negative charge Q_ as an approximation. In this case the integral 
Eq. can be evaluated imediately. If we further employ a spherical shell model for the 
nuclear charge distribution, a simple scaling-law is derived. This relates the energy shift 
Eq. ( p2p with the first-order Uehling correction AE\ p (Ueh) according to: 

AEj PVPb (^l,WK) ~ =— A£VP(Ueh) . (35) 

The validity of this scaling-law is supported by the fact that the numerical results for 
the Uehling correction do not depend significantly upon details of the extended nuclear 
charge distribution. With the scaling-law at hand we have an additional tool for testing the 
numerical results for the Wichmann-Kroll-in-Uehling contribution Eq. (|32| ) for strongly 
bound electrons. 

The complete evaluation of the Wichmann-Kroll-type contribution 

A £VPVPb ( ^ 2) = ( ^| W VPVPb |0A) (36) 

will be performed numerically according to the subtraction scheme introduced by Eq. 
([24]) respectively (p5|) . Fig. |] illustrates the subtraction scheme we applied for numerical 
calculations. 



VII. RESULTS AND DISCUSSION 



The scaling-law Eq. ( p5|) derived above for the "Wichmann-Kroll in Uehling" part 
AE\ pyph (J r l, WK) may only be useful as test for its complete numerical calculation. 
In order to employ the scaling-law the first-order Uehling correction ^^(Ueh) and the 
Wichmann-Kroll charge density p WK need to be determined. This can be achieved by 
means of very accurate numerical procedures developed earlier (see e.g. [13,15,3| and ref- 



erences therein). For a discussion of technical details encountered with the evaluation of 
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Pren we re f er to || . The total induced negative charge number Z^f^ is obtained by inte- 
grating the Wichmann-Kroll density from the origin r = up to r = r_ where changes 
its sign. Table | gives numerical results for the first-order Uehling correction, the induced 
negative charge number and for the corresponding energy shifts Ai?J PVPb (J r l, WK) ac- 
cording to the scaling-law for the lS-ground state of hydrogen-like Lead and Uranium 
respectively. The energy correction is repulsive, since the induced charge is negative. 
Let us compare the effect of the higher-order contribution to the two-loop correction 
Ai?J PVPb (h.o.) defined by Eq. (|30|) with the one obtained in Uehling approximation 
AE\ PVPh (J r l, Ueh). The latter is tabulated in the first column of Table for K- and L- 
shell electrons in hydrogen-like Lead and Uranium. This Uehling-in-Uehling contribution 
is attractive and amounts to less than 20% of the total Kallen-Sabry correction to the 
Lamb shift as tabulated in ||10|| . The results for the exact two-loop correction are tabulated 
in the second column of Table 

In comparison with the exact numerical results for the Wichmann-Kroll-in-Uehling con- 
tribution as presented in the third column of Table |TJ the scaling-law indeed gives the 
right order of magnitude but leads to a systematic overestimation due to the neglection 
of the long-range tail of The numerical results for the Wichmann-Kroll- type con- 

tribution A£j PVPb (T2), which have been obtained according to the representation (|25|), 
are presented in the last column of Table |I[ As for the Uehling-in-Uehling correction 
this contribution carries the same overall sign, i.e. it also acts attractive. It comprises 
about 40 % of the two-loop diagram evaluated in Uehling approximation. This result also 
indicates that higher-order (Za) contributions to the one-loop polarization insertion of 
the photon propagator may be not small compared with the Uehling approximation. Ac- 
cordingly, it may be not surprising at all, if the higher-order contribution to the S(VP)E 
- correction (see Fig. [I]) may turn out to be important as well. 

Even in the case of Uranium the order of magnitude for the total higher-order part 
A£ , v pvPb (h.o.) ( gee Table y) is about _ iq-2 e y for the is- s tate, which indicates that 

the effect of the higher-order contributions to the two-loop vacuum polarization alone are 
far too small to be detected via Lamb-shift measurements with current accuracies. Fur- 
thermore, this correction is about one order of magnitude below the natural limitations 



for tests of QED set by nuclear polarization effects |L6|. It might also be instructive to 
compare this effect with the uncertainties of the lS-energy level caused by the uncer- 
tainties in the determination of the nuclear radii. The rms-radius for Uranium is given 
by < r 2 >\j 2 = 5.8604(23) fm |17j and leads to an uncertainty of about 5E± S ~ 0.1 eV. 
Although a complete numerical evaluation of the exact Kallen-Sabry diagrams, i.e., the 
higher-order contributions of the self-energy-corrected one-loop diagram VPVPc remain 
to be performed, the goals of this paper have been achieved: 1. We have presented a 
calculational scheme for evaluating the two-loop vacuum polarization corrections to all 
orders in the interaction with the external Coulomb potential. The problem is reduced 
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to the evaluation of an effective one-loop correction. Thus, the partial wave decomposi- 
tion can be employed in a similar way as it has been used successfully in calculations of 
the ordinary first-order vacuum polarization correction. 2. Numerical results for exact 
two-loop vacuum polarization correction to the Lamb shift of K- and L-shell electrons in 
hydrogen-like Lead and Uranium have been obtained. It turns out to be below the natural 
limits set by nuclear polarization effects and by the uncertainties of nuclear parameters. 3. 
Although the effect of the higher-order contributions to the two-loop diagram turn out to 
be small a further uncertainty in Lamb-shift calculations has been eliminated. Aiming for 
a relative precision of 10~ 6 of theoretical predictions of the binding energy we conjecture 
that the still unknown exact two-photon self energy corrections remain the major source 
of uncertainties. 
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TABLES 

TABLE I. Estimate of the " Wichmann-Kroll in Uehling" contribution to the lS 1 / 2 -Lamb 
shift in hydrogen-like Lead and Uranium. The values given for the first-order Uehling correction 
i?^ p (Ueh) are calculated assuming a uniform sphere model for the nuclear charge density |t|]. 
The Wichmann-Kroll density is calculated for a spherical shell model §§. 



system 


A^P /2 (Ueh) [eV] 


^•WK 


Ai?VPVPb (HiWK) [eV ] 


238 tt 
92 u 

»Fb 


-93.58 
-50.70 


-0.006 
-0.004 


0.0061 
0.0024 
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TABLE II. The results for the two-loop vacuum polarization contribution in Uehling ap- 
proximation A^ PVPb (Jfl,Ueh) (Eq. (||)) in comparison with the exact two-loop correction 
A£^[ PVPb to the Lamb shift of strongly bound electrons in hydrogen-like Lead and Uranium. 
The higher-order contributions AEy™^, WK) (Eq. @) and A J BX PVPb (J r 2) (Eq. @) 
are listed separately. Energies are given in units of eV. 



system 


state 


A£X PVPb (^l,Ueh) 


A^PVPb 


AjB vpvPb^ 1)WK ) 


A#X PVPb (JR2) 




lsl/2 


-0.1150 


-0.1530 


0.0040 


-0.0420 


238 tt 
92 u 


2S1/2 


-0.0220 


-0.0286 


0.00074 


-0.0073 




2^1/2 


-0.0023 


-0.0036 


0.000079 


-0.0014 




lsi/2 


-0.0520 


-0.0685 


0.0015 


-0.0180 


20 8 8 2 Pb 


2«l/2 


-0.0092 


-0.0118 


0.00027 


-0.0029 




2pi/ 2 


-0.0006 


-0.0010 


0.000018 


-0.0004 
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FIGURES 
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FIG. 1. QED corrections of order a 2 in hydrogen-like atoms. The double lines indicate 
wave functions and propagators in the external Coulomb field of the nucleus. 
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FIG. 2. VP-dressed one-loop vacuum polarization potential (indicated by shadowed lines). 
The "R" inside of the external VP-loops refers to the renormalized one-loop VP-potential V^f 
induced by the external Coulomb field of the nucleus, which appears as additional external 
potential. 
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FIG. 3. F-dressed one-loop vacuum polarization potential (indicated by shadowed lines) 
according to the decomopsition (0). (y) symbolizes interactions with the total (renormalized) 
external potential Eq. (p^). 
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FIG. 4. Diagrammatic representation the subtraction scheme for deducing the renormalized 
two- loop vacuum polarization potential according to Eqs. ( p2| ) - (pip. The first term in curly 
brackets stands for the part Uj~f VPb and the second term in square brackets symbolizes the 
Wichmann-Kroll-type potential ll^pp^^- 
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FIG. 5. Diagrammatic representation the subtraction scheme for deducing the higher-order 
in (Za) part of the renormalized two-loop vacuum polarization potential U^^ Ph . The first term 
in curly brackets stands for the part U^Ywk an d the second term in square brackets represents 
the Wichmann-Kroll-type potential Z^|" VPb according to Eq. (^5|) . The third one-loop potential 
involves Dirac states in the presence of the vacuum polarization potential only. 
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FIG. 6. Uehling-vacuum polarization charge density ep^,„ h induced by the extended external 
charge of a Uranium nucleus (uniform sphere model) as a function of the radial distance r. 
Natural units are used. 
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